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Employing time-dependent projection formalism, a Fokker-Planck equation with non- 
Markovian transport coefficients is derived for large amplitude collective motion. Proper- 
ties of transport coefficients for diffusion processes in a potential well and along a potential 
barrier are discussed, and their connection with the fluctuation-dissipation theorem is in- 
vestigated. In the case of diffusion in a potential well, the formalism naturally leads to the 
generalization of the well known Einstein relation by including quantum fluctuations in 
addition to thermal fluctuations. Furthermore, it is shown that at low temperatures dif- 
fusion along potential barrier is significantly modified by quantum fluctuations. Explicit 
expressions for transport coefficients are presented for the case of the Caldeira-Leggett 
model. 



I. INTRODUCTION 

The fluctuation-dissipation theorem is a fundamental relationship in open system problems such as the 
Brownian motion (!]]. The most famous one is the Einstein relation, which connects the diffusion coefficient 
to the friction coefficient with the temperature as the proportional coefficient. This relationship can be 
easily derived if the bare potential for the Brownian particle, or more generally the collective motion, is 
fiat and if the environment can be approximated by a canonical distribution, ft has to be modified if the 
quantum fluctuation becomes comparable to the thermal fluctuation. When the bare potential for the 
collective motion is a parabolic well, the modified formula is well known. It resembles the original Einstein 
relation, in which temperature is replaced by an effective temperature pip!- The effective temperature 
converges to the genuine temperature when the thermal energy dominates the quantal zero point motion 
energy of the collective motion, and it converges to the zero point motion energy in the opposite limit. 

However, there exist only a few investigations when the bare potential of the collective motion is not a 
potential well but a potential barrier . The aim of this paper is to address to this question. Following 
the standard treatment for deriving the fluctuation-dissipation relation for motion in a potental well, 
we investigate the connection between moments of the response and the correlation functions for the 
motion along a potential barrier. We wish to emphasize that this study is important in various problems. 
One of the examples is the stability of nuclear collective excitation . Another example is the fusion 
between two very heavy nuclei, which is used to try to synthesize superheavy elements. In this case, the 
fission barrier locates inside the fusion barrier. Consequently, it is not sufficient to overcome the fusion 
barrier in order for the two nuclei to fuse. They have to overcome also the fission barrier. One possible 
theoretical approach to this problem is to apply the Fokker Planck equation, which has been originally 
put forward by Kramers || to describe nuclear fission and chemical reactions, to describing the time 
evolution from the fusion point to the region inside the fission barrier. Though several pioneering works 
in this method have already been published 0, it is not obvious whether one can apply the original 
Fokker Planck equation of Kramers as it is. For example, it is required to synthesize superheavy elements 
at temperatures as low as possible in order to sustain the shell correction energy which yields the barrier 
against fission by quantum effect. On the other hand, the curvature of the fission barrier is of the order 
of 1 MeV. One will then need to care about the quantum fluctuation. In this connection, it is highly 
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desirable to generalize the original Fokker Planck equation by Kramers by taking quantum fluctuation 
into accont for the diffusion process along a potential barrier. The formalism presented here provides a 
first step in this direction. Following a weak-coupling treatment, we derive a generalized Fokker-Planck 
equation which incorporates quantum statistical effects in terms of non-Markovian transport coefficients. 

The paper is organized as follows. In section II, we sketch the basic ideas of our theory to describe the 
time evolution of coupled collective and microscopic motions. It is a generalization of the quasi linear 
response theory in ref. ]8| to include a quantum effect of the collective motion. We thus derive the gen- 
eralized Fokker-Planck equation for the collective motion with non-Markovian transport coefficients. In 
section III, we discuss the properties of the force exerted by the microscopic motion, which we sometimes 
call the bath space, to determine the average classical trajectory of the collective motion. In section 
IV, we represent the response and the correlation functions in terms of the spectral functions ||, and 
discuss their symmetry properties. In section V, we discuss the memory times of the response and the 
correlation functions. In order to have some concrete ideas, we resort to the Caldeira-Leggett model jl0| 
to this end. Section VI deals with the collective motion in a flat potential or in a parabolic well and 
derives the well known fluctuation dissipation relation including the effect of quantum fluctuation. In 
sect. VII we consider the case of Ohmic dissipation of the Caldeira-Leggett model, and give the explicit 
expressions of the odd moment of the response function and the even moment of the correlation function 
for the motion in a potential well and the odd moment of the response function for the motion along a 
potential barrier. In sect. VIII we study the properties of the even moment of the correlation function 
for the collective motion along a parabolic potential barrier by taking the case of the Ohmic dissipation 
of the Caldeira-Leggett model with a cut-off function as an example. We show the strong time depen- 
dence of the even moment of the correlation function and show that the asymptotic value as the time 
goes to infinity converges to a classical expectation value at high temperatures, but significantly deviates 
from it at low temperatures. We summarize the paper in sect. IX. Appendix A gives the main steps 
of the perturbation treatment of the coupled von-Neumann equations for the collective and microscopic 
motions and show how the reduced von-Neumann equation with non-Markovian transport coefficients is 
obtained. Appendix B is added to briefly sketch the symmetry property of the spectral function which 
is useful to derive the fluctuation-dissipation theorem. Appendix C collects some basic formulae and the 
spectral function of the Caldeira-Legget model. In Appendix D we present a general expression of the 
m-th moment of the response function for the Caldeira-Leggett model. We give an explicit expression of 
the linear friction for the case of Ohmic dissipation. Also we show that the linear friction is absent in the 
case of the super-Ohmic dissipation, and that the leading order is the n-th oder friction if the spectral 
density is proportional to w™. 



We call the subspace of the collective motion and that of microscopic motions or environments as the 
subspaces A and B, respectively, and use A and B as the indices to distinguish the quantities referring to 
each subspace. The mass, coordinate and momentum for the collective motion are denoted by M, q and 

p, respectively, while those for the environments by m„ Xi,pi, i being 1,2 We often use x to denote the 

ensemble of the coordinates of the environmental degrees of freedom. As often done, quantum operators 
are distinguished from the corresponding c- numbers 0(t) with O. 

We assume that the Hamiltonian for the total system is given by 



II. BASIC EQUATIONS OF THE QUASI LINEAR RESPONSE THEORY 



A. Galilei transformation 



H = H A + H B + V c (q, x) 



(1) 



with 



f) 

H A = — + U(q) 
2M w 



(2) 
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where U(q) and V c (q,x) are the bare potential for the collective motion and the coupling Hamiltonian, 
respectively. The time evolution of the density operator of the total sysytem is determined by the von- 
Neumann equation, 



ih^- t p(t) = [H,p(t)} 



(3) 



We first introduce the classical variables q(t) and p(t) in the spirit of the Ehrenfest theorem by 

q(t)=Tr(qp(t)), p(t) =Tr(pp(t)). (4) 



Eq.|| leads to 



M ' 
1 



p(t) = -Tr(\p,H]p(t)) 
in 

= ^Tr({[p,U(q)] + [p,V(q,x)]}p(t)) 



(5) 



(6) 



We then introduce the Galilei transformation operator and the density operator in the Galilei transformed 
coordinate system by 



It is easy to prove that 



with 



G(t) = exp[-(p(t)q-q(t)p)} 
p{t) = G(t)p G {t)&{t) 



ihg- t PG(t) = [H G (t),p G (t)} 



(7) 
(8) 

(9) 



H G (t) = &(t)HG(t) - in&{t)-G{t) 



(10) 



The Galilei transformation is convenient to describe a large amplitude collective motion, because one 
then needs to handle only the fluctuations around the mean values. Noting that, 



&(t)qG(t)=q + q(t) 
G\t)pG{t)=p + p(t) 

we expand H G in the power series of the operators of the collective motion 

-ri>(<) • / \ i . r ./s dU(q(t)) dV c (q(t),x), 
H G = c- number +P{ P ±> - *(*)} + q{p(t) + + — g^l 



(11) 
(12) 



p- 1 d 2 U d 2 V n 
2M + 2^dq 2 + Ihf 



H B {x) + VM*),*) + ^T7 + + ^Ut)f + ^ (3) 



(13) 



where + ^r) q (t) is the value of + ^p 1 ) at ? = <?(£) and V^ 3 ^ is the third and higher order 
terms in fluctuating quantities. We note that by taking the time derivative of the equation 

Tr(p PG (t))=0 (14) 

we obtain 
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P(t) = -— -Tr( — p G (t))+T (15) 

where T is the force exerted by the third or higher order fluctuations and it reads 

1 r) 3 TI f) 3 V 

T= h Tr{{ w + W e ' G[t)}+ (16) 

We ignore the V*- 3 -* in eq.(p"3|) and correspon ding ly the force J- in cq.()l5|) in this work. On the other 
hand, the q 2 term on the right hand side of eq.(fL3[), which has been ignored in is essential to discuss 
the quantum fluctuation, which is the major issue in this paper. Using eqs.Q and ( fi"5"| ) the Hamiltonian 
in the Galilei transformed coordinate system can be written as 

v 2 1 

H G = c- number + q ■ F(q{t), x) + h B (t) + + -C(q(t))q 2 (17) 

where h B {t), the effective Hamiltonian for the subspace B at time t, the effective coupling force F(q(t), x) 
and the curvature parameter C(q(t)) are given by 

h B (t) = H B (x) + VMt),x) (18) 

F( q (t), x) = - Tr(«^p G (t)) (19) 

8 2 TI d 2 V 

C(l{t)) = ( W + -g^U) (20) 

The von-Neumann equation for pg{^) reads 

ih-p G (t) = [h B (t) + q ■ F(q(t),x) + ^ + -C(q(t))q 2 ,p G (t)]. (21) 



B. Fokker-Planck equation for the collective motion 

We now introduce the time evolution operator for the subspace B, u B , by 

ih^-u B (t,t ) = h B {t)u B (t,t Q ) (22) 
at 

with the initial condition 

u B (t ,t ) = l (23) 

and the density operator D(t) by 

PG (t) = u B (t,t )D(t)u B (t,t ) (24) 

The density operator D(t) obeys 



i^D(t) = [q ■ f(q(t),x),D(t)} + + \c{<lU) Uf . Din] (2:,) 



where 

f(t)=u B {t,t )Fu B {t,t ) (26) 
We introduce the reduced density operators for the subspaces A and B by 
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D A {t) =Tr B D(t), 
and the correlation density by 

D c (t) = D(t) - D A (t) ■ D B (t) 
These density operators obey the following coupled equations 



D B (t) = Tr A D(t) 



a 



ih^D A (t) = Tr B [q ■ f(q(t), x), D c (t)] + + -C(q(t))f , D A (t)] 
ih^t> B (t) = Tr A [q ■ f(q(t) ,x),D c (t)] 
ih^D c (t) = [q ■ f(q(t),x),D A (t) ■ D B (t)} + (1 - D A (t)Tr A - D B {t)Tr B )[q ■ f(q(t),x),D c (t)] 



(27) 
(28) 

(29) 
(30) 

(31) 



Following a similar treatment presented in ||, whose major steps are briefly described in the Appendix 
A, we derive a diffusion equation for the collective motion by handling eqs.(|2^) through (|l]) in the second 
order perturbation theory with respect to the fluctuating force /. In the derivation, we ignore the time 
dependence of the curvature parameter C, and introduce the moment functions C{t,t\) and S(t,t\) by 



-,(°), 



(t, ti)qu^(t, h) = C{t 7 ti)q - S(t, ti)p 



where 



They are 



U^(t,tx) = exp[^{^ + l -C{q{t))q 2 ){t ~ h)] 



C(t,ti) = cos[fl(t - ti)} 
S{t 7 tx) = 1 sin[n(t - t x ) 



with 



when C > and 



V M 



C(t,tx) = cosh[fl(t - tx)] 
S(t,tx) = , 1 sinhMt - tx) 

y/M\C\ 



with 



V M 

when C < 0. 

Furthermore we define the response and the correlation functions by 



(32) 

(33) 

(34) 
(35) 

(36) 

(37) 
(38) 

(39) 



X { J(t,tx) = ^Tr B ([f a (t),Mtx)]D B (tx)) 



X$(Mi) = ^Tr B ([f a (t),Mtx)] + D B (tx)) 



(40) 
(41) 
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and their moments by 



X% E) (*) = j dhC(t, h) X i-\t, h) (42) 
X { -°\t)= I dtiSfahfr^bh) (43) 

Jto 

X% E \t)= fdtiCfrtJxtyfrh) (44) 

Jto 

X%°\t)= fdtiS&h^frh) (45) 

Jto 

In these equations the lower indices a and (3 represent the spatial component in the case when the 
diffusion process takes place in a multi-dimensional space. The upper indices E and O stand for even 
and odd, respectively. The von-Neumann equation for £>A(t) is then given by 

[x% E \m f 3 - x%°\t)mD A {t)]] 

- ^[fe, lX ( ~f \t)q p - X t-°\t)pp,b A [t)} + ] 

(46) 

Taking the Wigner transform of each side of eq. ( J4q ) we obtain the Fokker Planck equation for the Wigner 
distribution function for the subspace of the collective motion Daw (<JS P, t) @ defined by 

/OO 11 

dre -W« < q + -r\D A (t)\q - -r > (47) 

The result reads 

+ <I7^» 0, 55s + ^sfe'^''" (48) 



This is a generalization of the well known Fokker Planck equation in literatures [plJllfl . Although it appears 
to be a Markovian equation, the memory effects are effectively incorporated in the transport coefficients. 
As a result, the formalism provides a basis for describing quantum statistical effects of collective motion. 
We also note that, since the derivation of eq.(^8|) is based on a perturbation theory with respect to the 
fluctuating force /, the formalism is valid in the weak coupling limit, where the time scales of the response 
and correlation functions, i.e. their memory times ri ' and ri + ' , are much shorter than the time scale, 



i.e. the relaxation time tr, of the collective motion fll2| , |13| . In sect.V, we discuss the properties of the 
response and correlation functions for the Caldeira-Leggett model, and show that their memory times 
Tc ' and Tc + ^ get very small if the cut-off frequency is large. This is consistent with the well known fact 
that the memory time gets very short if the collective motion couples with many incoherent degrees of 
freedom. 



C. Coupled equations to determine fluctuations of the collective motion 

We note that the average values of q and p are zero in the Galilei transformed space, and that the 
solution of eq.(|4q) is a Gaussian. Therefore, one can put 



G 



Daw(p, ?>*) 



(2tt) 3 / 2 A 1 /2 



- y 



2A 



where yi = q and 1/2 = P and the 2x2 matrix er is the inverse matrix of the 2x2 matrix 



(49) 
(50) 

(51) 



which determines the fluctuations, i.e. the mean square deviations from the average values. The values 
of <xy are obtained by solving the following coupled equations given by eq. (f48|) , 



it 1 M*> 

<7 pp (t) 







2_ 

M 



M 





M 



-2{C- X ^ E) ) 



M 



<J qq (t) 
<7pp(t) 





2x t + E) 



(52) 



D. The Wigner distribution function in the original space 

The Wigner distribution function for the collective motion in the original space fixed frame is given by 

PAw(q,P,t) = D AW (q - q(t),p- p(t),t) (53) 

once the Wigner distribution function in the Galilei transformed space Daw is obtained from eqs.(^9|) 
and (H). 

III. CLASSICAL TRAJECTORY OF THE COLLECTIVE MOTION 

The classical trajectory of the collective motion given by q(t) and p(t) are determined by eqs.(|5|) and 
(|l5|). Let us express the force induced by the coupling to microscopic motions in terms of the response 
function. It is useful, for example, to discuss the fluctuation-dissipation theorem. We first note that to 
the same order as eq.(fl6|) the density operator for the subspace B, /5b (t) = Tr A p(t), obeys the following 
von-Neumann equation 

ih^p B (t) = [h B (t),p B (t)] (54) 

We determine ps (t) by a perturbation theory up to the second order of the coupling V c . The effective 
force exerted on the collective motion from the bath space is then given by 

dV c (q(t),x) 

F B A = -Tr B { q PB\t)) (55) 

oq 

= -Tr B ( dVci f ^ e A**(«-fe)fe(to)e-**Bft-«°)) 
oq 

_ I (\ tlTrB{ ?YMh^ [e ^HB^^^ (56) 

J to vq 
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A. Limit of ignoring the quantal effect 

There are two alternative ways to write down the connection between the friction tensor and the 
moments of the response function. The standard way is to ignore the quantum or the memory effect and 
relate the friction coefficient to the first moment of the response function. As often done, in this approach, 
we approximate /3b (to) by a canonical distribution pg (t) = e^*)^ - ^ 5 ) , ignore the time variation of f3, 
and expand V c (q(ti), x) in the square bracket on the r.h.s. of eq.(p6|) around t\ — t. We then obtain 

F BA = -Tr B ( dV ^ £) pp\t)) 

i% J dq 
-q(t)x ( - 1) (t) + (57) 

where the first moment of the response function X (*) * s defined by 

X ( ~p\t)= I dti(*-ti)xS(*,*i) (58) 

J t 

Eq. ( |57j ) is the expansion with respect to the moments of the response function, whose higher order terms 
can be ignored if the memeory time is short enough, ft shows that the linear friction coefficient is given 
by the first moment of the response function in such cases. Note, however, that it becomes zero in some 
cases. As we show in Appendix D, the super-Ohmic dissipation is one of such examples. 



B. Inclusion of quantum effect 

The second way is to take the quantum or the memory effects into account and represent the friction 
coefficient in terms of the odd moment of the response function x^~°^ to be consistent with eq. (p|) . In 
this method, we first remark the following relationship between the classical variables at two different 
times, 

q(h) = C(t, t x )q{t) - MS(t, h)q(t). (59) 

This is valid when the classical motion evolves along a parabolic potential well or a parabolic barrier 
without friction ,i.e. in the harmonic approximation. We expand V c (q(ti),x) in the last term on the 
r.h.s. of cq.(|56|) in powers of q(h) — q(t) around q(t), and use eq. (|5^) . This is allowed to the lowest order 
of the coupling strength. We then obtain 

FBA = -Tr B ( dV ^ £) p^Ht)) 

ih J dq 
+ q(t) [ tfti[C(t,*i)-l]x (_) (*.*l) 

Jto 

-Mq{t)x { ~°\t) + (60) 

We note that the odd moment of the response function xl~a°' ) (*) converges to the first moment of the 
response function in the limit of small Q, i.e. in the limit of small C. More precisely, for small ft 

(-o)m- J W- x) (*)-o(n 2 ) (oo) 



1 lx { - l) {t) + o{n 2 ) (c<o) (bl) 



The results in eqs.(p7|) and (|60| ) agree to each other in this limit. 
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IV. PROPERTIES OF THE RESPONSE AND THE CORRELATION FUNCTIONS IN THE 

LINEAR RESPONSE THEORY 



A. Spectral function 

We first note that TrB(f a (t)fp(ti)DB(ti)), which appears in the definition of the response and the 
correlation functions, can be rewritten as 

TrB(f a (t)f (h)DB(h)) = Tr B {F a {t)u B {tM)h{t^B{tM)pB{t)) (62) 

where 

mb (Mi) = wb(Mo)«b(*i>*o) (63) 

Similarly, 

Tr B (Mti)f a (t)D B (ti)) = Tr B {uB{tM)h(ti)u B {tM)F a {t)pB{t)) (64) 

We consider the collective motion in the weak coupling limit, for which tJ ,+ ' << tr, so that the 
change of the collective variable during the decay time ri of the collision kernel is small. Then, wc 
can approximate these relations by 

Tr B {f a {t)ff}{ti)D B {ti)) « Tr B {F a {t)u B {tM)Mt)u B {tM)pB{t)) (65) 
TrB(U(ti)f a (t)D B (h)) « Tr B {uB{tM)h{^B{tM)F a {t)p B {t)) (66) 

where the time evolution operator Us(Mi) is given by 



u B 



m tt V h B (t')dt' 

(Mi) =Te sfi J *i 1 (67) 



with the time ordering operator T. For time intervals much shorter than the relaxation time of the 
collective motion, it can be apparoximated as, 

u B {t, ti) = u B (t; t = t - ti) w e w &B(t)T (68) 
We now define the spectral functions by 



(69) 
(70) 



/-t-00 
dTe WT Tr B {F a {t)u B {t; T)F (t)u B {t; r)p B (t)) 
-OO 

/ + °° 
dTe WT Tr B (F a (t)e-* hB ^ T Ff}{t)e* hB ^ T PB {t)) 
-OO 

Inversion of these equations yields, 

/+00 j 
— e- iuT J a ^(t;w) (71) 



OO 
+ OO 



Tr B {uB{t-T)Fp(t 1 )iJ B {t-T)F a {t)p B (t))=l Jp, a (t;u) (72) 
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B. Symmetry properties 

We introduce a further approximation by replacing the density operator of the subspace B with a 
canonical distribution, 

p B {t) ~ exp{j3(t){F-h B (t))} (73) 
As we show in Appendix B, the resultant specral function has the following symmetry property 

J , a (t;u,)=e^ hu J a , l3 (f,-u;) (74) 

One can then rewrite eq.(|7^) as 



/ + OG J 
™-i»r e -lHt)*»j pfaw) (75) 
-oo 2lT 

Consequently, the response and the correlation functions can be expressed as 

X { -l{t,h)=^-kt-,T = t-H)= / — e -^-{l-e-^}J a , p (t^) (77) 



Also, one can prove the following symmetry properties 

x$(*;-t) = x$(*;t) (78) 

X i a :i(t;-r)=~x ( a :l(t;r) (79) 

Furthermore, if we define 



xi + l (*,«>)= dr cosut x« (t ; r) (80) 
Jo 

poo 

xkl( t l u; )= dr sinujT x ( al{t]T) (81) 



then we obtain 



X { +l(t;u) = - J a , a (t\u) (82) 

X ( aMt;u) = ±r{l-e-W R "}J a<a (t;u) (83) 



V. MEMORY TIMES OF THE RESPONSE AND CORRELATION FUNCTIONS 

An important issue to further proceed is the memory times of the response and the correlation func- 
tions. In order to see their connection to various quantities of the system such as the properties of 
the environments, or the coupling Hamiltonian and the temperature, we consider explicitly the case of 
Feynman- Vernon model JlJ], which is a popular model in open system problems and has been widely 
used to discuss the macroscopic quantum tunneling since the seminal work of Caldeira and Leggett |l(| . 
Here we discuss the case of Ohmic dissipation. It is easy to show that the response function is given by 

2 df f°° 

Xi:l(t;r) = -(£) 2 r, sin(wT)Gjg(uj) du> (84) 
ir dq J 
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where g(uj) is the cut off function (see Appendix C). Using its explicit form, we have 



( _j ^ | 1 ^(iL) 2 ri ■ ojI ■ t ■ e~i w ^ 2 (Gaussian cut-off) 
1 (^) 2 'J' li, c 'e _WcT (Drude regularization) 

Eq.(^||) shows that the memory time gets shorter for larger cut-off frequency. On the other hand, the 
correlation function is given by 



X&fo r) = ^(|) 2 -T-f du, cos{ur) ■ £coth( 



P(t)hw 



) • g(u) (86) 



Figure 1 shows the correlation function X (t> t) as a function of r for three values of the cut-off 
frequency and for a given temperature T=l MeV. Since we consider a linear coupling, X (*i T ) i s 
independent of time t, so that we denoted the ordinate as X ( T ) instead of X (^i T )- The correlation 
function is normalized to 1.0 at r = 0. It decreases rapidly with r and converges to zero after passing 
through a minimum as r approaches infinity. Figure 2 shows the correlation function as a function of r 
for three values of temperature and for huj c = 20.0 MeV . It exhibits a similar behavior as in figure 1. 
We define the memory time rj as the time interval during which the correlation function becomes e _1 
of its value at r = 0. The horizontal dashed line in Figs.l and 2 denotes the position of e _1 . Figure 3 
shows the deduced memory time of the correlation function as a function of the cut-off frequency for three 
temperatures studied in Fig. 2. The memory times for T=0.1 MeV and 1.0 MeV are indistinguishable 
in the figure. Figure 4 shows the memory time as a function of the temperature for three values of the 
cut-off frequency. These results show that the memory time is nearly independent of temperature T and 
the cut-off frequency u c provided Tiuj c > 5 MeV, and its magnitude is around ri ~ (0.5 — 3) * 10~ 22 
second. Consequently, when the temperature is comparable or even lower than the energy of the quantum 
zero point motion, a sufficiently large cut-off frequency guarantees a short memory time even at such low 
temperatures. 



VI. FLUCTUATION-DISSIPATION THEOREM FOR A FLAT POTENTIAL OR A 
POTENTIAL WELL IN THE SHORT MEMORY TIME APPROXIMATION 

We are now ready to derive the fluctuation-dissipation theorem in the case, where C > 0, i.e. for the 
diffusion process in a flat potential or in a potential well. We first change the integration variable in 
eq.(^) from t\ to r = t — t% 

xL 0) (*) = / ° drS(t; t) X U (*5 r) (87) 
where the quantity S(t; r) is given by, 



o 



S(t, ti) = S(t; T = t-t 1 )= 1 sm[Or] 

v MC 

We assume that the memory time of the response function ri is much shorter than the relaxation time 
of the collective motion tr. This condition is satisfied in the weak coupling limit if the collective motion 
couples with many incoherent environmental degrees of freedom as mentioned before, or if the cut off 
frequency is high in the case of Caldeira-Leggett model as shown explicitly by eq.(^). Then, replacing 
the upper limit of the integration in eq.([37|) by oo, the odd moment of the response function becomes, 



1 r°° 

xL 0, (t)" f dTsinSlTxtt(t,T) 



(89) 



_^_L {1 _«,-«<)-} ,„(«;„) (91) 
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Assuming that the memory time of the correlation function is also much shorter than the relaxation 
time of the collective motion, i.e. by assuming Tc + ^ << tr, we obtain 

i i i p -p(t)m 

X { a + a E Ht) = Ja,a(t;n) (92) 

Eqs.(|9l|) and d92| ) lead to the well known quantum fluctuation-Dissipation relation, 

Y (+E) (t) 

*^W=MT* (93) 

with the effective temperature given by 

T* = -Ml ■ coth{-j3(t)hQ) (94) 



2 v 2' 

t (t » nn) 
\nn (t « mi) 



(95) 



VII. DISSIPATION COEFFICIENT FOR THE CALDEIRA-LEGGETT MODEL WITH A 

CUT-OFF FREQUENCY 

In this section, we derive explicit expressions of the odd moment of the response function for the case 
of Ohmic dissipation of the Caldeira-Leggett model. 



A. Collective motion in a flat potential or in a potential well 

We first consider a collective motion in a flat potential or in a potential well. Using the definition of 
the odd moment of the response function given by eq.([43"|) and the expressions derived in sect.V for the 
response function, we obtain for the Ohmic disspation 

xL 0) (t) = j;(f q ) 2 -V9m (96) 

Note that the magnitude of the cut-off function g(£i) is almost one as long as the cut-off frequency is 
much larger than the curvature parameter fl of the bare potential of the collective motion. On the other 
hand, the corresponding Xaa (t) is given by 

xi + a E) (t) = V ■ (|) 2 • ■ • coth{\l3{t)m (97) 



Eqs.(p6|) and fl97|) guarantee the fluctuation-dissipation theorem eq.(93). 



B. Collective motion along a potential barrier 

We next consider a collective motion along a potential barrier, where the curvature parameter C is 
negative. In this case, eqs.([43|), ( |38|) and fl84| ) lead to 

2 1 df /" t_t o r°° 
xL 0) (t) = -V / — - (^) 2 / drsinhiVLT) dwsimuT ■ u ■ g(w) (98) 
7T y/M\C\ aq J J 

for the case of Ohmic dissipation. For the Gaussian cut-off and the Drude reguralization, it becomes, 
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f -L^ffWf' ta drsinhinr) . T . e -\^^f (Gaussian cut-off) 

xi~fHt)=l ^ m Ct _ (99) 

| n— 7==(^) 2 ti)^ J Q drsinh^flr) ■ e (Drude regularization) 

Furthermore, by replacing the upper limit of the r integral with oo, we obtain 

(-0)/ t j = I -V-^ )2 2 (Gausszan cut-off) 

I 3f(^) 2 "f I^nfe+oy (£>™cte regularization) 

where w c > has been assumed for the case of the Drude regularization. These expressions resemble the 
result for the collective motion in a flat potential or in a potential well given by eq.(|9^). It is, however, 
interesting to notice that the effect of the frequency cut-ff appears in a slightly different way. 

The results in this section seem to indicate that the dissipation coefficient is almost independent of the 
temperature of the environmental space. However, this should not be taken as a general conclusion, but 
rather as a feature of the linear harmonic oscillator coupling model. For example, the friction coefficient, 
which is calculated in a random matrix model in p"5| for describing energy dissipation in deep inelastic 
heavy-ion collisions, rapidly increases with temperature. 



VIII. CONNECTION BETWEEN THE MOMENTS OF THE RESPONSE AND 
CORRELATION FUNCTIONS FOR THE DIFFUSION ALONG A POTENTIAL BARRIER 



Similar expression to ec 
eqs.dH), ©,© and " 



can be obtained for the moment of the correlation function by using 



It reads 

h df ft— to r<x> -i I p -f3{t)hw 

xi + a £) « = -^) 2 y o drco8h(Qr)J du coswt ■ U • g(u) ■ j^^- (101) 

h df o 1 f^* f°° 1 + e ~ fJ ^ huj 

= 2 J dT cosh ^ J ^ duJ C0SUJT ' w ' ' 1 _ e -/?(t)fto, ( 102 ) 

n ,df i r*-* u(c\ \ f°° a ( ^ +hS S^Lx f in ,i 

= — rj( — ) — / dr coshlilT) / dui cosuit ■ ui ■ g(ui) ■ coth{ ) (10o) 

7T dq 2J J_ oc 2 

We first perform the r integral analytically then the u) integral using Cauchy's residue theorem. The 
time integral introduces a pole at u> = ±«0, and the last factor introduces series of poles associated with 
the Matsubara frequency (2tt /Tij3) ■ ni, ni being an integer. There appears another pole at u> — ±iui c if 
we use the Drude regularization for the high frequency cut-off. The resultant formula reads 

xi^M— ^) 2 (n + r 2 + y 3 ) (104) 



with 



7r dq ' 



Y x = --u? c ■ e — «(*-*o) • cot (-p{t)nu Jc ){—^—e nl - t - t ^ + -J— e -^-t o)} (105) 
4 2 uj c — \l uj r + \i 



Y 



2 



\ p fi/3(t) v l L °> . £ 



' e n ^ 1 ^ + e -"(*-*o)} (106) 



7r-2n o " 7r-2n , o 



Yl = ''~^J^'r l ' r,ll: ^ j{i ( 11,71 



Note that the Yi term should be ignored in the limit of large u> c . If lo c » T » ft, then keeping only 
n=l term we have 
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(108) 



and 



Y * ~ m (109) 

If we keep only the I3 term in this case, asymptotically, i.e. as time goes to infinity, we have 

*str>«— )~4 ,2 w (110) 

which together with eq.( |100| ) leads to the classical fluctuation-dissipation theorem. 

The even moment of the correlation function (t) 1S the momentum diffusion coefficient. At 



low temperatures, a significant deviation in the diffusion coefficient from the classical value eq.(llC) is 
expected from quantal effects ||. In order to see this explicitly, in figure 5 we illustrate the normalized 
diffusion coefficient Y(t) — [Yi(t) + Y%{t) + 13(f)] — as a function of time f for several temperatures for 
the collective frequency Ml = 1.0 MeV and the cut-off frequency Tvjj c — 20.0 MeV. The horizontal lines 
show the classical values determined by temperature Y c i as = T. At high temperatures, the normalized 
diffusion coefficient monotonically increases with time and gradually approaches to its classical value. 
However, the behavior of the diffusion coefficient is very different at low temperatures. After a sharp rise, 
it undershoots the classical value and may even become negative for large times. We note that a similar 
behaviour of the diffusion coefficient has been discussed also in ||] in connection with non-Markovian 
effects. 



IX. SUMMARY 



We discussed the quantum effects on the transport coefficients for a large amplitude collective motion. 
To that end, we derived a Fokker Planck equation with non-Markovian transport coefficients based on 
a time dependent perturbation theory. We studied both cases, where the collective motion moves in a 
flat potentail or in a potential well, and along a potential barrier. We have shown that our formalism 
naturally leads to the known generalized fluctuation-dissipation theorem which combines the thermal and 
quantal fluctuations when the bare potential is a potential well. In the case when the collective motion 
moves along a potential barrier, the well known classical behaviour of the moment of the correlation 
function which governs the fluctuation is realized only at high temperatues. Using a numerical evaluation 
of the moment of the correlation function, we demonstrated that it significantly deviates from the classical 
expectation at low temperatures, both with respect to the time dependence and the asymptotic value. 
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APPENDIX A: DERIVATION OF THE REDUCED VON NEUMANN EQUATION IN THE 
GALILEI TRANSFORMED COORDINATE SYSTEM 

In this Appendix, we explain the major steps to derive eq. (|46|) from eqs.(|2^) through (EH]) by the second 
order perturbation theory. We first multiply a parameter A to the fluctuating force /, and use it as a 
parameter of the perturbation theory. 
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As usual in the perturbation theory, we expand the reduced density operators as 

oo 

D A (f) = 5> n £<r>(t) (Ai) 

n=0 

oo 

D B (t) = J2^ i B ) (t) (A2) 

oo 

Dc(t) = J2^c l) (t) (A3) 

n=0 

To the zeroth order with respect to A, we have 

in^\t) = [^+ 1 -C(t)f,vf(t)} (A4) 



whose solutions are 



^ 0) (t) = #i 0) (t)6i%)[^°W (A7) 
V^(t) = T>^(t ) (A8) 
P^(t) = C/f (t)^° 3 (t )[C/W(t)]t = (A9) 



U A °\t)=T{exp[^ f {f- p + \c{t x )q 2 )dt l ]} (All)) 



with 



L ih J to 2M 2 

where T is the time ordering operator. We have assumed that the two subspaces A and B are decoupled 
at t = t , so that T>£\to) = 0. 

To the 1st order with respect to A, we have 

ih^-V^it) = Tr B [q ■ f(t),V§\t)] +[^+ \c{t)?M\t)) (All) 
^'(t) = Tr A [q ■ mV$\t)] (A12) 
ih^\t) = [q ■ f{t),f>f{t)f>f{t)] + [q ■ f(t),V<§\t)] 

-T)f{t)Tr A [q ■ /(t), - vf{t)Tr B [q ■ f(t),V$\t)} 

+[^ + \c{t)e,V^\t)}, (A13) 



where cq.( A13) can be solved by the method of variational constant by putting 

V%\t) = uf\t)C{t)\uf\t)y (A14) 

The result reads 

+ ±rf dt 1 {U A a \t,t 1 )q{U%\t,t 1 tf ■ /(iO.^tMOPfliij^lMi)]^^)] (A15) 
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itli 



Now we note that we have 



i^-Dfit) = Tr B [q ■ f(t),D^(t)} + [JL + ic(t)«f , 



in the second order approximation. Note that 



t)f{t) = f>f{t) + XD<?(t) + x 2 vf{t) 



i(o) 



Eq.( [A15| ) then leads to 



ih^T)f(t) = i / * 1 Tr s [^/(t),[C/i 0) (t,t 1 )g[i/7(t,t 1 )]t./(t 1 ), 
/to 



r<°>/ 



"9t ih ' A 



1 



dtiTr B [? • /(f), [^(MiMt/rtMi)] 1 • D^W^Cti)]] 



7t(°) 



n(2)wn(2), 



to 



(A16) 



(A17) 



(A18) 
(A19) 



(A20) 



(A21) 



The transformation from eq.( A20|) to eq.(A21) is consistent with the second order perturbation theory. 
It is now straightforward to transform eq.(A21) into the form of eq. (Eg) . 



APPENDIX B: SYMMETRY PROPERTY OF THE SPECTRAL FUNCTION 



The spectral function is given by 

f+OO 



(Bl) 



for a slow collective motion. We assume that the density operator of the bath space p B (f) can be 
approximated by a canonical distribution e^ t ^~ hB ^ t '> . The spectral function then reads 



+00 



J a ^(t;u) = / dTe toT Tr B (^ a (t)e-* AB W r F /9 (t)e* fcfl W T e /3 W^- fca( * )) ) 



The key step is to insert two delta functions as follows 



— 00 
00 



J 0l/ j(f;«)= / dTe WT Tr B {F a {t)e-* hB ^ T dES(E - h B {t))F {t) 

> J — OO 



We then obtain 



/+00 r+00 
dE j dE'6(E - E' - huj)j a43 {E, E') Peq (E') 
-OO J — CO 



(B2) 



(B3) 



(B4) 



1G 



where 

j a AE,E') =Tr B (F a (t)8(E-h B (t))F p (t)8(E'-h B (t)) (B5) 
Peq (E)=e^- E *> (B6) 

One can easily prove the following symmetry relations 

j (3ia (E',E)=j atl3 (E,E') (B7) 
\j atl3 {E,E')}* =j F t <F t(E',E) (B8) 

We now consider 

/+OO /- + 00 

dE / dE'S{E - E' - huj)j fj , a {E, E') Peq {E') (B9) 
-oo J — CO 



Using eqs.(B7) and (B6) and interchanging the integration variables E and one can prove that 

Jf,, a (t; u) = e mn "J a At; ( B1 °) 
This proves eq.([74|) in the main text. 

APPENDIX C: BASIC FORMULAE FOR THE FEYNMAN-VERNON-CALDEIRA-LEGGETT 
MODEL IN THE LINEAR RESPONSE THEORY 

In this Appendix, we collect several basic formulae for the case, where the environment is represented 
by an aggregate of harmonic oscillators, and where the collective motion couples them linearly with 
respect to their coordinates 10 The unperturbed Hamiltonian of the bath space H B , the coupling 



Hamiltonian V c (q(t),x and the force operator F a (q(t),x) read 

HB=J2(^Pi+l m ^ ( C1 ) 

i 

V c (q(t),x) = f(q(t))J2ciXi (C2) 

i 

F a (q(t),x) = ^M)ly2c i {x i -Tr(x i p B (t))} (C3) 
dq a 

In the linear response theory, we replace the p B (t) and h B (t) on the r.h.s. of eq.([70|) by the canonical 
distribution e 0{t){F-H B ) anc j jj B , respectively. Eqs.(|Cl]) through (|C3|) then lead to 

fjf ft p -/3(t)Huii i 

4° a LLR \t; U ) = 2, ■ {j- q f ■ £ c?— + u*) ^—^ + S(u - ^——Q-} (C4) 

One can easily prove that this expression satisfies the symmetry property eq.(|74|). It would be physically 
reasonable to assume that all u>i are positive. Then, for to > 0, the relevant part of the spectral function 
is 

^->( ti „) = 2 ..(|) 2 .2 c? ^ J („-^_i^ <„>„) (C5) 

We now define the spectral density J(ui) by 

J( u ) =y^(% S(u-Ui) (C6) 

. rmuJi 
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Following 0, we assume that 



J(w) - -r)—^ I u n g(u) (C7) 

7T Wo 

Th e ca ses of n=l and n=3 correspond to the Ohmic and super Ohmic dissipations, respectively. In 
eq.(C7),we have introduced ujq in oredr to keep the dimension of 77 independent of n. g(tu) is the cut-off 
function, for which we assume two-types of functional forms 

g{u) = gi(u) = e- ( - )2 (C8) 
9 {uj) = 92 {uj) = 1 (C9) 

We call the first and the second the Gaussian cut off and the Drude regularization, respectively. 



APPENDIX D: GENERAL EXPRESSION OF THE M-TH MOMENT OF THE RESPONSE 
FUNCTION FOR THE CALDEIRA-LEGGETT MODEL 



We consider 



k aq u jo 



dTT T 



dujsiriDT ■ uj n ■ g(tu) 



(Dl) 



where both n and m are integers; n,m=l,2,.... This will be the m-th moment of the response function if 
the spectral density is given by 



(D2) 



One can show 



Epecially, 



K dq oj, 
' dq x % 



duj 1 



l(l) 2 ^T(-l) lfl+1 (-l)' 



x[- 1 \t)=v(f q ) 2 9(0) 



"duo 1 



-{co n g(u)}]^ Q 



(D3) 
(D4) 

(D5) 
(D6) 



Eq.(D5) agrees with the results eq.(61) together with eqs. (|96|) and (100). On the other hand, eq.(D6) 
means that there exists no linear friction in the case of super Ohmic dissipation. Eq.(D4) shows that the 
lowest order friction is n-th order friction if the spectral density is proportional to oj n as given in eq.(D2) 
and the friction coefficient is given by the n-th moment of the response function. 
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Figure captions 

Figure 1. The normalized correlation function x^ + H T ) as a function of r for three values of the cut-off 
frequency and for a temperature T=l MeV. The horizontal dashed line shows the value e _1 . 

Figure 2. The normalized correlation function x^ + H T ) as a function of r for three values of temperature 
and for a cut-off frequency hui c — 20.0 MeV. The horizontal dashed line shows the value e" 1 . 

Figure 3. Memory time of the correlation function ri + ' as a function of cut-off frequency Tiu) c for three 
values of temperature. 

Figure 4. Memory time of the correlation function ri + ^ as a function of temperature for three values of 
the cut-off frequency Huj c . 

Figure 5. The normalized diffusion coefficent Y(t) as a function of t for several values of temperature 
and for a collective frequency Ml = 1.0 MeV and a cut-off frequency Uuj c = 20.0 MeV. The horizontal 
dashed lines show the corresponding classical values. 
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